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Abstract 

We modify the usual Erdos-Renyi random graph evolution by letting connected 
clusters 'burn down' (i.e. fall apart to disconnected single sites) due to a Poisson 
flow of lightnings. In a range of the intensity of rate of lightnings the system sticks 
to a permanent critical state. 



Keywords: forest fire model, Erdos-Renyi random graph, Smoluchowski coagulation 
equations, self- organized criticality 

AMS subject classification: Primary 60K35 

Acknowledgement. This research was partially supported by the OTKA (Hungarian 
National Research Fund) grants K 60708 and TS 49835. The authors thank the anony- 
mous referee for thoroughly reading the manuscript. Her/his comments and suggestions 
helped us improving the presentation. 

Submitted to EJP on September 1, 2008, final version accepted May 14, 2009. 



1 



Contents 



Introduction 

1.1 Context 

1.2 The model 

1.3 The main results 

Coagulation and fragmentation 

2.1 Forest fire flows 

2.2 The Markov process 

2.3 The integrated Burgers control problem . . . 

Boundary behavior 

3.1 Elementary facts about generating functions 

3.2 Bounds on E 

3.3 No giant component in the limit 

The critical equation 

4.1 Elementary properties 

4.2 Bounds on E' 

4.3 Uniqueness 



1 Introduction 
1.1 Context 

In conventional models of equilibrium statistical physics, such as Bernoulli percolation, 
random cluster models, the Ising model or the Heisenberg model there is always a pa- 
rameter which controls the character of the equilibrium Gibbs measure: in percolation 
and random cluster- type models this is the density of open sites/edges, in the Ising or 
Heisenberg models the inverse temperature. Typically the following happens: tuning the 
control parameter at a particular value (the critical density or the critical inverse temper- 
ature) the system exhibits critical behavior in the thermodynamical limit, manifesting 
e.g. in power law rather than exponential decay of the upper tail of the distribution of 
the size of connected clusters. Off this particular critical value of the control parameter 
these distributions decay exponentially. We emphasize here that the critical behavior is 
observed only at this particular critical value of the control parameter. 

As opposed to this, in some dynamically deflned models of interacting microscopic 
units one expects the following robust manifestation of criticality: In some systems dy- 
namics defined naturally in terms of local interactions some effects can propagate in- 
stantaneously through macroscopic distances in the system. This behavior may have 
dramatic effects on the global behavior, driving the system to a permanent critical state. 
The point is that without tuning finely some parameter of the interaction the dynamics 
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drives the system to criticality. This kind of behavior is called self-organized criticality 
(SOC) in the physics literature. The two best known examples are the sandpile mod- 
els where so called avalanches spread over macroscopic distances instantaneously, and 
the forest fire models where beside the Poissonian flow of switching sites/edges from 
"empty" to "occupied" state (i.e. trees being grown), at some instants connected clusters 
of occupied sites/edges (forests of trees) are turned from "occupied" to "empty" state 
instantaneously (i.e. forests hit by lightnings are burnt down on a much faster time scale 
than the growth of trees). These models and these phenomena prove to be difficult to 
analyze mathematically rigorously due to the following two facts: (1) There are always 
two competing components of the dynamics (in the forest fire models: growing trees 
and burning down forests) causing lack of any kind of monotonicity of the models. (2) 
Long range effects due to instantaneous propagation of short range interactions are very 
difficult to be controlled. 

Regarding forest fire models there are very few mathematically rigorous results de- 
scribing SOC. The best known and most studied model of forest fires is the so-called 
Drossel-Schwabl model. For the original formulation see [TT], or the more recent survey 
[T6] . We formulate here a related variant. 

Let A„ := 1/ n [—n,nY. The state space of the model of size n is Vtn '■= {0, 1}^": 
sites of An can be occupied by a tree (1) or empty (0). The dynamics consists of two 
competing mechanisms: 

(A) Empty (0) sites turn occupied (1) with rate one, independently of whatever else 
happens in the system. 

(B) Sites get hit by "lightnings" with rate X{n), independently of whatever else happens 
in the system. When site is hit by lightning its whole connected cluster of occupied sites 
turns instantaneously from "occupied" (1) to "empty" (0) state. (That is: when a tree is 
hit by lightning the whole forest to which it belongs burns down instantaneously.) 

The dynamics goes on indefinitely. 

As long as n is kept fixed the mechanism A + B defines a decent finite state Markov 
process - though a rather complicated one. The main question is: what happens in the 
thermodynamic limit, when n oo, A„ y Z'^? Can one specify a dynamics on the 
state space floo '■= {0, 1}^ which could be identified with the infinite volume limit of the 
systems defined above? 

In order to make some guesses, one has first to specify the lightning rate A(n). In- 
tuitively one expects four regimes of the rate A(n) with essentially different asymptotic 
behavior of the system in the limit of infinite volume: 

I. If A(n) <^ |A„|~^ then the effect of lightning is simply not felt in the thermodynamic 
limit: in macroscopic time intervals of any fixed length no lightning will hit the entire 
system. Thus, in this regime the system will simply be the dynamical formulation 
of Bernoulli percolation. 

II. If A(n) = |A„|~^A with some fixed A G (0, oo) then one expects in the thermody- 
namic limit the following dynamics (described in plain, non-technical terms). The 
system evolves as dynamical site percolation, with independent Poisson evolutions 
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on sites, and with rate \0(t), where 6(t) is the density of the (unique) infinite clus- 
ter, the sites of this (unique) infinite cluster are turned from occupied to empty. 
After this forest fire the system keeps on evolving like dynamical percolation until 
a new infinite component is born, and the dynamics goes on indefinitely. 

III. If |A„|^^ ^ \{n) <^ 1 then in the infinite volume limit - if it makes any sense - 
something really interesting must happen: The lightning rate is too small to hit finite 
clusters within any finite horizon. But it is too large to let the infinite percolating 
cluster to be born. One can expect (somewhat naively) that in this regime in the 
thermodynamic limit a dynamics will be defined on ^loo in which in plain words the 
following happens: 

- empty (0) sites turn occupied (1) with rate one, independently of whatever else 
happens in the system; 

- when the incipient infinite percolating cluster is about to be born, it is switched 
from "occupied" (1) to "empty" (0) state; 

- the dynamics goes on indefinitely. 

In this way this presumed infinitely extended dynamics would stick to a permanent 
critical state when the infinite incipient critical cluster is always about to be born, 
but not let to grow beyond criticality. 

IV. If \{n) = A G (0, oo) then lightning will hit regularly even small clusters and thus, 
one may expect that - if the infinitely extended dynamics is well defined - the system 
will stay subcritical indefinitely. 

There is no problem with the mathematically rigorous definition of the infinitely 
extended dynamics in regimes I. and II. But these plain descriptions don't necessarily 
make mathematical sense and it is not at all clear that such infinitely extended critical 
forest fire models can at all be defined in a mathematically satisfactory way. 

In our understanding, the most interesting open questions are the existence and char- 
acterization of the infinitely extended dynamics in regime III. and/or the A — oo limit 
in regime II. and/or the A ^ limit in regime IV., after the thermodynamic limit. 

There are however some deep results regarding these (or some other related) models 
of forest fires, though clarification of the above questions seems to be far out of reach at 
present. 

Here follows a (necessarily incomplete) list of some important results related to these 
questions: 

- M. Diirre proves existence of infinitely extended forest fire dynamics in a related model 
in the subcritical regime IV. , [l2]. In a companion paper he also proves that under 
some regularity conditions assumed the dynamics is uniquely defined, [13]. 

- J. van den Berg and R. Brouwer, respectively R. Brouwer consider the so called self- 
destructive percolation model, which is very closely related to what we called regime 
II. above. They prove various deep technical results and formulate some intriguing 
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conjectures related to the A — > oo limit in regime II. (of the already infinitely extended 
dynamics), see [2], [3], [8] 



- J. van den Berg and A. Jarai analyze the A ^ asymptotics of the (infinitely extended) 
model in regime IV. in dimension 1, [Ij. 

- J. van den Berg and B. Toth consider an inhomogeneous one dimensional model which 
indeed exhibits SOC, see |5]. (In one dimensional space-homogeneous models of course 
there is no critical behavior) 



We investigate a modification of the dynamical formulation of the Erdos-Renyi ran- 
dom graph model, adding "forest fires" caused by "lightning" to the conventional Erdos- 
Renyi coagulation mechanism. Actually our model will be a particular coagulation- 
fragmentation dynamics exhibiting robust self-organized criticality. 

Let Sn '■= {1,2, ... ,n} and i3„ := = ■ i,j G iS„, i ^ j} be the set of 

vertices, respectively, unoriented edges of the complete graph /C„. We define a dynamical 
random graph model as follows. The state space of our Markov process is {0, l}'^". 

Edges of /C„ will be called occupied or empty according whether uj{i,j) = 1 or 
uj{i,i) = 0. As usual, we call clusters the maximal subsets connected by occupied edges. 

Assume that initially, at time t = 0, all edges are empty. The dynamics consists of 
the following 

(A) Empty edges turn occupied with rate 1/n, independently of whatever else happens 
in the system. 

(B) Sites of JCn get hit by lightnings with rate A(n), independently of whatever else 
happens in the system. When a site is hit by lightning, all edges which belong to its 
connected occupied cluster turn instantaneously empty. 

In this way a random graph dynamics is defined. The coagulation mechanism (A) 
alone defines the well understood Erdos-Renyi random graph model. For basic facts and 
refined details of the Erdos-Renyi random graph problem see [H], [6], [I5]. As we shall 
see soon, adding the fragmentation mechanism (B) may cause essential changes in the 
behavior of the system. 

We are interested of course in the asymptotic behavior of the system when n ^ oo. 
In order to formulate our problem first have to introduce the proper spaces on which our 
processes are defined. 



1.2 The model 



We denote 



V: = {v= (i;,)^^^ : Vk>0, J]t;fc<l}, ^(v) := 1 - J] t;,, 



(1) 




Vi :={v G V : ^(v) = O}. 



(2) 
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We endow V with the (weak) topology of component-wise convergence. We may interpret 
as the density of the giant component. 

A map [0, oo) 3 v(t) G V which is component- wise of bounded variation on com- 
pact intervals of time and continuous from the left in [0, cxd), will be called a forest fire 
evolution (FFE). If v(t) G Vi for all t G [0, oo) we call the FFE conservative. Denote the 
space of FFE-s and conservative FFE-s by £, respectively, £i. The space £ is endowed 
with the topology of component-wise weak convergence of the signed measures corre- 
sponding to the functions Vk{-) on compact intervals of time. This topology is metrizable 
and the space £ endowed with this topology is complete and separable. 

Now, we define the cluster size distribution in our random graph process as follows 

Vn,k{t) '■= E Sn '■ j belongs to a cluster of size k at time t} =: rr^Vn,k{t)-, (3) 

v„(t) := K,fc(t)),gi,. (4) 

This means that v„(t) is the cluster size distribution of a uniformly selected site from 
Snt at time t. Clearly, the random trajectory t i— > v„(t) is a (conservative) FFE. We 
consider the left-continuous version of t v„(t) instead of the traditional c.a.d.l.a.g., for 
technical reasons discussed in Subsection 12. 1[ 

We investigate the asymptotics of this process, as n — >■ oo. 

It is well known (see e.g. [9], [lO], [1]) that in the Erdos-Renyi case - that is: if 
\{n) = 

Vn(-) v(-) = (t;fc(-))fcgN n^oo, (5) 

where the deterministic functions t i— > Vk{t) are solutions of the infinite system of ODE-s 

k^ 

^fc(^) = 2 ^ vi{t)vk-i{t) - kvk{t), k>l, (6) 
1=1 

with initial conditions 

Vk{0) = 4,1. (7) 

The infinite system of ODE-s (l6|) are the Smoluchowski coagulation equations, the initial 
conditions ((71) are usually called monodisperse. The system ([6]) is actually not very scary: 
it can be solved one-by-one for = 1, 2, . . . in turn. For the initial conditions ([7]) the 
solution is known explicitly: 

Vk{t) = ^e-'H'-\ 

(^fe(^))fcLi G V is a (possibly defected) probability distribution called the Borel distri- 
bution: in a Galton- Watson branching process with offspring distribution POI{t) the 
resulting random tree has k vertices with probability ffc(t). Thus the branching process 
is subcritical, critical and supercritical for t < 1, t = 1 and t > 1, respectively. 
For general initial conditions ffc(O) satisfying 

oo oo 

Y^VkiO) = l, ^A;V(0)<oo, 

k=l k=l 
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the qualitative behavior of the solution of ([6j) is similar: Define the gelation time 

oo 

Tgei:= (8) 
fc=i 

- For < t < Tgei the system is subcritical: 6'(v(t)) = and, k i-^ ffc(t) decay exponen- 
tially with k. 

- For Tgei < t < oo the system is supercritical: 9{'v{t)) > and k Vk{t) decay 
exponentially with k. Further on: t \—>- 0(v(t)) is smooth and strictly increasing with 
limt_oo^(v(t)) = 1. 

- Finally, at t = Tgei the system is critical: 6'(v(Tgei)) = and 

oo 

J^t;/(rgei) X k^oo. (9) 

l=k 

Our aim is to understand in similar terms the asymptotic behavior of the system 
when, beside the Erdos-Renyi coagulation mechanism, the fragmentation due to forest 
fires also take place. 

Similarly to the Drossel-Schwabl case presented in subsection [Hi] we have four regimes 
of the lightning rate A(n), in which the asymptotic behavior is different: 

I. : X{n)<t:n~\ 

II. : X{n)=n-^X, Ag(0,oo), 

III. : n-^ < X{n) < 1, 

IV. : A(n) = A G (0,oo). 

The n — » oo asymptotics of the processes t ^— v„()f:) in the four regimes is summarized 
as follows: 

I. The effect of lightnings is simply not felt in the n — >^ oo limit. In this regime the 
system will be the dynamical formulation of the Erdos-Renyi random graph model, 
the asymptotic description presented in the previous paragraph is valid. 

II. In the — > oo limit the sequence of processes t i— > v„(t) converges weakly (in 
distribution) in the topology of the space £^ to a process t v(t) described as 
follows: The process t v(t) evolves deterministically, driven by the Smoluchovski 
equations (l6|) (exactly as in the limit of the dynamical Erdos-Renyi model) with the 
following Markovian random jumps added to the dynamics: 

P ( v(t + dt) = Jv I v(t) = V ) = Xe{\)dt + o{dt) (10) 
where J : V ^ V, (Jv)^ = t;^ + 5fc,i^(v). (11) 
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In plain words: with rate A6'(v(t)) the amount of mass 6'(v(t)) contained in the gel 
(i.e. the unique giant component) is instantaneously pushed into the singletons. 

III. This is the most interesting regime and technically the content of the present paper. 
In the n — i> oo limit (JS]) holds, where now the deterministic functions t Vk{t) are 
solutions of the infinite system of constrained ODE-s 

k ''^^ 

'"kit) = IT 5Z ^'(^)^fe-i W - kvkit), k>2, (12) 



1=1 



J2vk{t) = l, (13) 

km 

with the initial conditions ((71). Mind the difference between the system (EI) at one 
hand and the constrained system (fT2ll + (fT3l) at the other: the first equation from ((61) 
is replaced by the global constraint (fT3l) . A first consequence is that it is no more 
true that the ODE-s in (fT2l ) can be solved for A; = 1,2,..., one- by-one, in turn. The 
system of ODE-s is genuinely infinite. Up to Tgei the solutions of (JE]), respectively, 
of (fT2|l + (fT3l) coincide, of course. But dramatic differences arise beyond this critical 
time. We prove that the system (fT2l) + (fT3l) admits a unique solution and for t > T„ei 



TT 
l=k 



J2vi{t)r^J^k-'/', as k^oo, (14) 



where [Tgei, oo) Bin-* ip{t) is strictly positive, bounded and Lipschitz continuous. 
This shows that in this regime the random graph dynamics exhibits indeed self- 
organized critical behavior: beyond the critical time Tgei it stays critical for ever. 
The unique stationary solution of the system (fT2l) + (fT3l) is easily found 



/2n - 2\ 1 „ 1,^/9 
V n - 1 / n V47r 

IV. In the n — * oo limit ([5]) holds again, where now the deterministic functions t i— > Vk{t) 
are solutions of the infinite system of ODE-s 

k °° 
^'fcW = 9Zl^'W^fc-iW-fc^feW-Afct;fc(t) + A4,i^/t;Kt), k>l, (16) 



1=1 



with the initial conditions in Vi. The system (|T6l) is again a genuine infinite system 
(it can't be solved one-by-one for /c = 1, 2, . . . in turn). The Cauchy problem (fTBl) 
with initial condition in Vi has a unique solution, which stays subcritical, i.e. for 
any t G (0, oo) k ^— Vk{t) decays exponentially. The unique stationary solution is 
closely related to that of ( fTSll : 

A2 ^ ' 



Vx,kioo) = (A + 1) ( 1 - (^^^^y ] ^fc(oo) 
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1.3 The main results 

We present the results formulated and proved only for the regime III: <^ A(n) <^ 1, 
which shows self- organized critical asymptotic behaviour. The methods developed along 
the proofs are sufficient to prove the asymptotic behaviour in the other regimes, described 
in items I, II and IV but we omit these (in our opinion less interesting) details. 

Theorem 1. If the initial condition \{0) G Vi is such that J2T=i ^^"^^(0) < +oo, andTg^i 
is defined by (I8|) then the critical forest fire equations (fT2l) + (fT3l) have a unique solution 
with the following properties: 

1. For t < Tgf.1 the solution coincides with that of 

2. For t > Tggi there exists a positive, locally Lips chitz- continuous function Lp such that 

viit) = -vi{t) + ip{t) (17) 

and (fT4l) holds. 

Theorem 2. Let P„ denote the law of the random FFE of the forest fire Markov chain 
v„(t) with initial condition v„(0) and lightning rate parameter ^ ^{^) ^1- If 
Vn(0) — >■ v(0) e Vi component-wise where Y^'^Li^^'^kit) < -\-oo then the sequence of 
probability measures P„ converges weakly to the Dirac measure concentrated on the unique 
solution of the critical forest fire equations ffT2l) + (fT3l) with initial condition v(0). In 
particular 

V£ > 0, t > lim P( \vn,k{t) -Vk{t)\ > = 

2 Coagulation and fragmentation 
2.1 Forest fire fiows 

In this section we investigate the underlying structure of forest fire evolutions arising 
from the coagulation-fragmentation dynamics of our model on n vertices. 

We define auxiliary objects called forest fire fiows: let qn,k,iit) denote times the 
number of {k, Z)-coagulation events (a component of size k merges with a component of 
size /) up to time t. Let Tn^kif) denote ■ k times the number of fc-burning events (a 
component of size k burns) up to time t. For the precise definitions see (|271) . (l28l) . (l30l) 
and fl3T|l . 

In Subsection 12.11 and Subsection 12.21 we precisely formulate and prove lemmas based 
on the following heuristic ideas: 

• The state v„(t) of the forest fire process on n vertices (see ([4])) can be recovered if 
we know the initial state v„(0), and the fiow: qn,k,iit) for all fc, / G N and rn^kif) for 
all k. The precise formula is (fTOll . 



9 



• (fTOll is similar to the equations (fT6l) . This will help us proving Theorem O if 1 <^ n 
and ^ \{n) <^ 1 then the random forest fire evolution Vn(t) "almost" satisfies 
the equations ffT2l) + (fT3l) that uniquely determine the deterministic limiting object 
v(t). We essentially prove that (fT2l ) is satisfied in the n — ^ cxo limit in Proposition 
[U of Subsection I2.2[ We prove that (fT3l) is satisfied in the limit in Subsection I3.3[ 



We define the moments of v G V as 

oo oo oo oo 

mo = ^Vk, mi = ^k-Vk, m2 = ^k'^-Vk, m3 = ^k^-Vk 

k=l k=l k=l k=l 

By (H]) and dSD mo = 1 if and only if v G Vi. 

Fix T G (0, oo). A map [0,T] 3 t v(t) G V is a a forest fire evolution (FFE) on 
[0, T] ifffc(-), k G N is of bounded variation and continuous from the left in (0,T]. Denote 
the space of FFE-s on [0, T] by £^[0,T] and the space of FFE-s with initial condition 
v(0) = V G V on [0, T] by £^v[0,T]. Note that a priori 6{-) = 1 — X]fceN'^'i-(') ii^ed not be 
of bounded variation. 

If v„(-) G £^[0, T] is a sequence of FFE-s then we say that v„(-) — > v(-) if Vn,k{,-) =^ Vk{-) 
for all G N where denotes weak convergence of the finite signed measures on 
[0, T] corresponding to the functions fn,fc(-) and ffc(-). Note that we did not require the 
convergence of 6'„(-) to 9{-). 

This topology is metrizable and the spaces £^[0,T] and £^v[0,T] endowed with this 
topology are separable and complete (by Fatou's lemma, lim„^ooV„(t) stays in V). 

Denote N := {1, 2, . . . } and N := N U {oo}. 

A forest fire flow (FFF) is a collection of maps [0,r] 3 t ^^ (q(t),r(t)) where for 
< s < t < T 

= gfc,z(0) < qk,i{s) < qk,i{t)^ = {^IkAt)) k,ieW ^'^''(^) = ^nAt), 

= rfc(O) < rfc(s) < rfc(t), r(t) = (rfe(t))fceN, n{t) = 

We define 

qk{t) ■.= J2lkAt), lit) ■= X]?fc(t), r{t) := Y^rkit) (18) 

im ken ken 

and assume the finiteness conditions q{T) < +oo, r(T) < +oo. All functions involved 
are continuous from the left in (0,T]. This is why we have chosen to consider the 
left-continuous versions of these functions rather than the traditional c.a.d.l.a.g.: the 
supremum of increasing left-continuous functions is itself left-continuous, thus the left- 
continuity of gfc, 9 and r automatically follows from the left-continuity of qk,i and r^. 

We say that the FFF [0, T] 9 t i— > (q(t),r(t)) is consistent with the initial condition 
v(0) = V G V if t v(t) defined by 

k 

Vkit) =Vk{Q) + -Y^ qi,k-i{t) - kqkit) - Tkit) + l{fe=i}r(t), ken. (19) 

(=1 
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is in ^v[0,r]. That is: for all t G [0,T] and A; G N Vk{t) > and EfceN^fc(^) ^ 1 holds. 
In this case we say that the FFF (q(-),r(-)) generates the FFE v(-). 

We denote by J-'v[0, T] the space of FFF-s consistent with the initial condition v(0) = 
V G V. For any v G V, JFv[0,T] 7^ 0, since the zero flow is consistent with any initial 
condition. 

At this point we mention that later we are going to obtain a FFF (q„(-), r„(-)) 
from a realization of our model on n vertices by (|27I) . (l28l) . (l30l) and (l3TI) . There is a 
FFF corresponding to the limit object as well: for the solution of the critical forest flre 
equations (fT2ll + (fT3ll (the uniqueness of the solution is stated in Theorem [T]) we define 
(q(-),r(-)) by 

qk,i(t) = Vk(t)vi(t), qoo,k{t) = qoo,oo{t) = 0, rfc(t) = 0, roo{t) = '^{t) (20) 

with the ifit) of (fTTll . In Definition [T] we define a topology on the space of FFFs. In later 
sections we are going to prove that 

(qn(-),r„(-)) (q(-),r(-)) 

from which Theorem [2] will follow. 

Summing (fTOll for A; G N we obtain a formula for the evolution of 0{-) defined in ([T]): 
for s < t 

K 00 

9{t) = 9{s) + lim ^ ^ ■ - QkA^)) + 

k=l l=K-^k+l 

00 

^ k ■ (gfc,oo(t) - qk,oc{s)) - (roo(t) - roo(s)) (21) 

k=l 

Later we will see that the term limi^^oo J2k=i ^'i^K-k+i ^ " {qk,i(t) — qk,i{s)) does not 
vanish for the FFF defined by (EDD for the unique solution v(t) of (pD + iini) if Tgei < s < t: 
this phenomenon is a sign of self-organized criticality. 

If (q(-),r(-)) is a FFF then the functions qk^i, qk, q, i^k and r (where k,l & N) are 
continuous from the left and increasing with initial condition 0: such functions are the 
distribution functions of nonnegative measures on [0, T]. By g(T) < +00 and r(T) < +00 
these measures are finite. We denote by the weak convergence of measures on [0, T], 
which can alternatively be defined by point-wise convergence of the distribution functions 
at the continuity points of the limiting function. 

Definition 1. Let (q„(-), r„(-)) = {{qn,k,i{-))k,ien Arn,k{-))km) ' n = 1,2,... be a se- 
quence of FFFs. Define qn,k{-), qn{-) and rn{-) for all n by (fTHll . 
We say that (q„(-), r„(-)) (q(-),r(-)) as n ^ 00 if 

wk,ieN qn,kA-) ^ qkA-) 

V A; G N g„,fc(-) qk{-) 

qn{-) q{-) 

WkeN r„,fc(-) ^ rfc(-) 
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Note that we do not require r„ oo(') =^ ^oo(') and gn,fc,oo(') =^ <lk,oo for k & N. Nev- 
ertheless these "missing" ingredients of the limit flow (q(-),r(-)) of convergent flows are 
uniquely determined by the convergent ones if we rearrange the relations ( fTSll : 



9fc,oo(^) := Qkit) - ^ qk,ii^), (22) 
/eN 

r-ooW :=r(t) -^rfc(t), (23) 

fceN 

goo,oo(t) := q{t) -2j2(lk{t) + J2 ^'^.'(^)- (24) 

ken fc,/eN 

In fact, r„^oo(-) 7^ '"oo(-) and qn,k,Qo{-) 7^ Q'fc,oo have a physical meaning in the forest flre 
model if (q„(-), r„(-)) is deflned by ( 1301) and (pB : 

• In the X{n) = 0{n^^) regime = g„,fc,oo(-) 7^ Q'fc,oo(-) ^ indicates the presence of 
a giant component. The precise formulation of this fact for the Erdos-Renyi model 
is (1361). 



If X{n) <^ 1 then only "large" components burn. Indeed in Proposition [T] we are 
going to prove that for all G N r„ fc(-) converges to in probability an n ^ oo. 
Thus by (l23ll we have r(-) = roo(-) in the limit. But Theorem [T], Theorem [2] and 
((201) imply that = r„,oo(i) ^ r^{t) = ip{s) ds > for t > T^^i- 



jFv[0,T] endowed with the topology of Definition [U is a complete separable metric 
space: 

Lemma 1. // (q„(-), r„(-)) e ^v[0,T] for all n e N and (qj-),rn(-)) ^ (q(-),r(-)), 
then (q(-),r(-)) e^v[0,T]. 

Proof. By the definition of weak convergence, qk,i, qk, q, rk, r are increasing left-continuous 
functions with initial value 0. We need to check that the functions Too, g'fc.oo, and goo.oo 
(defined by ([23l) . (l22l) and (l24l) . respectively) are increasing. We may assume that < 
s <t <T are continuity points of g^, '"fc and r for all fc, Z G N. 
By Fatou's lemma we get 

roo{t) - roo{s) = lim (r„(t) - r„(s)) - lim (r„,fc(t) - r„,fc(s)) 

n— >oo ' ' n—>-oo 

fceN 



> lim sup r„(t) - r„(s) - ^ (r„,fc(t) - r„,fc(s)) 
'^^^ V fceN 

= lim sup (r„,oo(^) - rn,oo{s)) > 0. 
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One can prove similarly that qk^oo is increasing for A; G N. In order to prove that 
qoo,oo{t) - goo,oo(s) > limsup (g 

n, 00,00 (t) - gn,oo,oo(s)) 

n^oo 

let an,k,i ■= qn,k,i{t) - qn,k,i{s) for k,l eN. By ([24]) we only need to check 

lim an,k,i - limsupa„,oo,oo > 2 ^ lim ^ a„,fc,/ - lim an,k,i- (25) 

n— >00 ' ' r)— >nn n~-*QO ' ' ' ' n^oo 

k,ien km leN k,i£N 

Let 

Km '■= {{k, I) : {k>m and I = m) or {I > m and k = m)} U {(m, 00)} U {(00, m)}. 

The left hand side of (EHD is lim inf„^oo ZlmeN f^n,m, the right hand side is Y^men li^^n^oo Pn,m, 
where Pn,m '■= J2{ki)eKm ^n,k,i, and the inequality follows from Fatou's lemma. 

Now that we have proved that the limit of convergent flows is itself a flow, we only 
need to check that the limit flow is consistent with the initial condition v, but this follows 
from the facts that £'v[0,T] is a closed metric space and the mapping from jFv[0,T] to 
£^v[0,T] deflned by (fTOll is continuous with respect to the corresponding topologies. 

□ 

Finally we deflne the space of all FFF-s as follows: 

D[0,T]:={(v,q(-),r(-)) : v G V, (q(-), r(-)) G ^"^[0, T]}. 

This space is again a complete and separable metric space if we define (v„, q„(-), r„(-)) 
(v, q(-),r(-)) by requiring v„ ^ v (coordinate-wise) and (q„(-), r„(-)) (q(-),r(-)). 

Lemma 2. For any C < 00 the subset 

ICc[0,T] :={(v,q(-),r(-))GP[0,r] : q{T) < C} 

is compact in VlO, T]. 

Proof. 

r K K 
_ k=l 1=1 k=l 

by q{T) < C, dominated convergence and qk,i = qi,k- Thus summing the equations ([I9] 
with coefficients | we get 

k 



lim 

K^oo 



1 1 



k=l k=l k=2 

The inequalities 



r(T)<2 + C, r^(T)<l + lc, r,(T) < (1 + |)^ (26) 

follow from v(T) G V and q(T) < C . 

By Kelly's selection theorem and a diagonal argument we can choose a convergent 
subsequence from any sequence of elements of /Cc'[0,T] with the limiting FFF itself being 
an element of /CcfO, T]. 

□ 
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2.2 The Markov process 

It is easy to see that in order to prove Theorem [2] we do not need to know anything 
about the graph structure of the connected components: by the mean field property of 
the dynamics the stochastic process v„(t) defined by ([3]) and ([4]) is itself a Markov chain. 
The state space of the Markov chain t V„(t) is: 

:= {V = {Vk)keN ■■ Vk e {0, k, 2k,...}, ^ l^fc = n] 

k>l 

The allowed jumps of the Markov chain are described by the following jump trans- 
formations for i < j: 

aij : {V G fi„ : Vi{Vj - jl{i=j}) > O} ^ 

(o-ijV)^ := Vk - il{fc=i} - il{fc=i} + {i + j)l{k=i+j}, 

Ti-. {V eQn-Vi> 0} ^ Qn, (riV)^ := Vk + il{k=i} - «l{fe=i} 
The corresponding jump rates are an,i,j, i>n,i '■ K_|_: 

a„«(V) := ((l + l{,=,})n)~V,(F,-jl|,=,}), 6„,,(V) := X{n)V. 
The infinitesimal generator of the chain is : 

Lnfiy) = $^an«(V)(/K,V) - /(V)) + ^6„,(V)(/(r,V) - /(V)). 

i<j i 

We denote by Qn,k,i(t) and by Rn,k(t) the number of 0"^ ^ -jumps, respectively fc-times 
the number of -jumps occurred in the time interval [0, t]: 

QnAiit) ■= (1 + Mk=i}) ■\{se[0,t] : V„(s + 0) = {ak,iVn) - 0)} | , (27) 
Rn,k{t) ■.= l{kj^i}k-\{s e[0,t] : V„(s + 0) = (rfcV„)(s-0)}| . (28) 



Finally, the scaled objects are 





''K,fe(t), 




v„(t) : = 


^ {Vn,k{t))^^^, 


(29) 




^Qn,k,lif)^ 


<ln,k,oo{t) = 0, 




■ {(ln,kAt)) k,leN' 


(30) 


rn,k{t) := 


^Rn,k(t), 


'"n,oo(^) = 0) 


rn(t) : = 


i''n,kit))km 


(31) 


Now, given T G (0, 


oo) and some 


initial conditions 


Vn(0) = V„ 


G Vi, clearly t h- 


v„(t) G 



Vi is a conservative FFE, generated by the FFF (v„, q„(-), r„(-)) G V[0,T] through (fTOll . 
We denote by P„ the probability distribution of this process on X'fOjT]. We will always 
assume that the initial conditions converge, as n — ^ cxd, to a deterministic element of Vi: 

lim Vn,kiO) = Vk, V := {vk)keN e Vi. (32) 
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Proposition 1. The sequence of probability measuresWn is tight onV[0,T]. If X{n) ^ 1, 
then any weak limit point P of the sequence P„ is concentrated on that subset ofV[0,T] 
for which the following hold for k, I G N; 



qk,ii^) = / Vk{s)vi{s)ds, qk(t) = / Vk{s)ds, q{t) < t, rk{t) = (33) 
Jo Jo 

v(0) = V. (34) 



Proof. There is nothing to prove about the initial condition (IMll : it was a priori assumed 
in d^. 

In order to prove the validity of the integral equations (l33l ). note first that it is 
straightforward that the processes qn,kA^), {qn,k,i){t), qn,k(t), {qn,k)(t), rn,k{t), (r„,fc)(t), 
defined below are martingales: 



qn,k,l{t) '■= qn,k,l{t) - I Vn,k{s)Vn,l{s)ds -\ i^Al v^ i^[s)ds, 

Jo ^ ^0 

~ ^9 l|j!j-,^R -|- 2]l|^=n 2,klLff^—i\ 

{qn,k,i)(t) := qn,k,i{'t) / Vn,k{s)vn,i{s)ds ^ / Vn,k{s)ds, 

n Jo n Jo 

qn,ki^) ■= qn,ki^) - / Vn,kis)ds+- / Vn,kis)ds, 



n 



_ I r 2k I 

i,k){t) := g„,fc(t)^ / {Vn,ki-'^y + Vn,kis))ds + — Vn,k{s)ds, 

n Jq n - 

1 

qn{t) ■.= qn{t) - t + - / mn,i{s)ds, 

n In 



{qn){t) := qn{tf - - U + I '^Vn,k{sYds\ + ^ I mn,i{s)ds 
^ \ "'o fc=i J n Jq 



Jo 

{rn,k){t) ■■= rn,k{tf - I Vn,k{s)ds. 

n In 
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From Doob's maximal inequality it readily follows that for any k,l eN and e > 

lim P( sup I g„,fc,/(t) - / Vn,k{s)vn,i{s)ds | > £:) = 0, 

n->oo 0<t<T Jo 

lim P( sup I qn,k{t) - / Vn,k{s)ds \ > e) = 0, 

0<t<T Jo 

lim P( sup qn{t) — t > e) = 0, 

n^oo 0<t<T 

lim P( sup I rnk{t) I > = 0. 
Hence (l33ll . Tightness follows from 

E(g„(T))<T, (35) 

Markov's inequality and Lemma [2l □ 

If we consider the case X{n) = (this is the dynamical Erdos-Renyi model) then 
(I5]) + (16]) follows from Proposition [T] since (fTOl l becomes 



Vk{t) = Vk{0) + -y2qi,k-iit) - kqk{t) = VkiO) + / -}^vi{s)vk-i{s) - kvkjs) ds 
^ 1=1 ^ 1=1 

which is the integral form of ([6]). Plugging ( 1331 ) into ( !22ll we get for t > Tgei 

gfc,oo(t) = f Vk{s)e{s)ds > 0. (36) 



2.3 The integrated Burgers control problem 

If v(-) e S^^[0,T] is generated by a FFF satisfying ((33|) through i^, then 

oo oo 

K-) = $^^;t(-) =5Z^fc(-)+^oo(-) =5Z0 + roo(-) =roo(-) 

fceN k=l k=l 

and v(-) is a solution of the controlled Smoluchowski integral equations with control 
function r(-): 



t^fe(t) = ^fc(0) + - V / i;Ks)^fc-/(s)rfs-A; / t;fc(s)rfs + l{fc=i|r(t), /c e N (37) 
^ Jo Jo 

oo 

t^fc(^)>0, 5^t^fc(t)<l (38) 
fc=i 

v(0) = vo G Vi. (39) 
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By g(T) < T, r^{-) = r(-) and ((26]) we get 

T 

= r(0) < r(s) < r(t) for < s < t < T, r(T) < 1 + -. (40) 

Using induction on k one can see that the initial condition vq and the control function 
r(-) determines the solution of ([371) . (l39l) uniquely. 
For V e V we introduce the generating function 

oo 

V : [0, oo) ^ [-1, 0], V{x) := J] ^;fce-^'^ - 1. (41) 

fc=i 

a; I— >^ V(x) is analytic on (0, oo) and has the following straightforward properties: 

lim V{x) = -1, V'{x) < 0, V"{x) > 0. (42) 

x— >oo 

It is easy to see that if 1 1-^ \(t) is a solution of ( l37l ). ( l38l ). ( l39ll then the corresponding 
generating functions t i— > V^(t, ■) will solve the integrated Burgers control problem 

V{t,x) -V{0,x) + [ V{s,x)V'{s,x)ds = e-''r{t), (43) 
Jo 

- 1 < ^(t, 0) < (44) 

V{0,x) = Vo{x). (45) 

The control function r(-) was defined to be continuous from the left in (fTSl) . but it need 
not be continuous: when A(n) = n~^X then the FFE obtained as the n ^ oo limit 
satisfies (1371 ). (1381 ). (1391) . but the control function r(-) evolves randomly according to the 
rules (IIOD, ([IH): 

P ( r(t + dt) =r{t)+ e{t) I T{t) ) = A^(t)rft + o{dt) 

Thus r(-) is a random step function in this case. 

In order to rewrite (Il3l) as a differential equation we introduce a new time variable r: 



t(r) :=max{t : t + r(t) < r} (46) 
It is easily seen that t(r) is increasing and Lipschitz-continuous: 

t(r) = /" a{s)ds <«(■)< 1 (47) 

Given a solution V{t,x) of ((43]), ((ll), (SSD define 

V(r, x) := F(t(r), x) + (r - t(r) - r(t(r))) e"" (48) 
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Then by S43h we have 

Y{t,x) = V{0,x) - V{s,x)V'{s,x)ds + {T-t{T))e-''. (49) 
Jo 

Now we show that for all r > 0, x > and t > we have 

drV{T, x) = -V(r, x) V'(r, x)«(r) + (1 - a(r))e-^ (50) 

-l<V(r,0)<0 (51) 

V{0,x) = Vo{x) (52) 

V{t + r{t),x) = V{t,x) (53) 

First note that the fact 

Y{t,x) ^V{t{T),x) =^ a(r)=0 (54) 

follows directly from ([46]), ((47]) and (SHI): if r(t+) ^ r(t), then a(r) = for all t + r{t) < 
T < t + r{t+). The differential equation ^ follows from ([47]), ^ and ((541). The 
boundary inequality (ISTl) follows from 



-1 < ■^^(t(r),x) < V(r,x) < \/(t(r) + ,x) < 0. 



The initial conditions (1451 ) and ( 1521 ) are equivalent, and (1531 ) follows from (1481 ) and (1461 ). 

From the definition of Lebesgue-Stieltjes integration it follows that for all ti < t2 we 
have t {t ) t 

/(t(r))(l - a(r)) dr = T f{t) dr{t) (55) 

'h+r{ti) Jti 

3 Boundary behavior 

3.1 Elementary facts about generating functions 

In this subsection we collect some elementary facts about generating functions, which 
will be used along the proof of Theorem [1] and Theorem [21 For v G V we introduce the 
generating function V{x) defined in ([411) which has the straightforward properties listed 
in (Il2l) . It is also easy to see that for any v G V and any x > 

\V'{x)\<^x-\ V^"(x)<0^x-^ \V"'{x)\<(^^^ x-^ (56) 

We define the functions E : (0, oo) — > (0, oo), E* : [0, oo) (0, oo], E^ : [0, oo) 
[0, cx)) as follows: 

^(^) ^= --^T^TTT' E*{x):= sup E{y), E^{x) := mi E{y) (57) 

V"[X) 0<y<x 0<y<x 

Note that these functions are continuous on their domain of definition. 



18 



Lemma 3. Let v G Vi. 

1. For any x > 

0< V{x)V'{x) <E*{x). (58) 

2. If in addition 

y'(0) := limF'(x) = -cx) (59) 
then the following bounds hold 

2V2^,(x)VV/2 < -Vix) <2'/^E*{xY/^x'/^ (60) 

2-'/^E,{x)E*{x)-'^^x~'/^ < -V'{x) < 2-^/^E*{x)E,{xY^/^x-^^/^ (61) 

2-^'^E,{xfE*{x)-'''^x-^/^ < V"{x) < 2-^'^E*{xfE,{x)-''l^x~'"^ 

E,{x) < V{x)V'{x) < E*{x). (62) 

Proof. Since v G Vi we have ^(0) = 0. Denote the inverse function of —V{x) by X{u): 
X{-V{x)) = X. Note that 

^'^) = j^^^y 

and thus 

X(0) = 0, X'(0) = -V'iO)-\ X'\u) = E{X{u))-\ 
It follows that for u G [0, -V{x)]: 

-V'{0)-^ + E*{x)-\ < X'{u) < -V'{0)-^ + E4xy\, 

-V\0)-'u + E%x)-\ < X{u) < -V\0)-'u + E,{x)-'-. 
Hence, all the bounds of the Lemma follow directly. □ 

3.2 Bounds on E 

We assume given a solution of the integrated Burgers control problem: ( i43ll . ( !44l ). (l45l l 
with a control function r(-) satisfying ( l40ll . 

We fix t G (0, oo), X G (0, oo). All estimates will be valid uniformly in the domain 
(t,x) G [0,t| X [0,x]. The various constants appearing in the forthcoming estimates will 
depend only on the initial conditions V{0,x) and on the choice of (t, x). The notation 

A{t,x) X B{t,x) 

means that there exists a constant 1 < C < oo which depends only on the initial 
conditions (l45l) and the choice of (t, x), such that for any {t,x) G [0,t| x [0,a;] 

C-^B{t, x) < A{t, x) < CB{t, x). (64) 
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The notation A(t,x) = 0{B(t,x)) means that the upper bound of (l6il) holds. 

In the sequel we denote the derivative of functions f(t,x) with respect to the time 
and space variables by f{t,x) and f'(t,x), respectively. 

First we define the characteristics given a solution of f l43l) . (i45l ). ( l44l ): for t > 0, x > 
let [0,t] 3 s ^t,x{s) be the unique solution of the integral equation 

^t,x{s) =x-V{t,x){t-s)+ [ (M-s)e-«*'^('')rfr(M). (65) 



Existence and uniqueness of the solution of (16511 follow from a simple fixed point argu- 
ment. Now we prove that (given (t,x) fixed) s C,t,x{s) is also solution of the initial 
value problem 

^6,x(s) =: itA-^) = Vis, 6,x(s)), 6,x(t) = X. (66) 

In order to prove this we define V(r, x) by fllSll . Thus from flSUl it follows that that the 
solution of (l66l l satisfies 

^6,.(t(r)) = V^(t(r), 6,.(t(r)))a(r) = V(r, 6,.(t(r)))«(r) (67) 



From this and ( 15011 we get that 

-^V(r,6,.(t(r))) = V(r,6,.(t(r)))+V'(r,6,.(t(r)))-;^6,.(tW) = (1 - a(r))e-«'-(*M) 
Integrating this and using ^t,x{t) = x and ( l53l l we get for all ri < t + r(t) 

pt+r{t) 

V(n,et,..(t(ri))) = K(t,x) - / (1 - a(r))e-f*-W^» rfr 



n 



Substituting this into the r.h.s. of (l67ll . integrating and using (mil we get for all T2 < 

t + r{t) 

itM^2)) = x- V{t, x){t - tir,)) + / (t(r) - t(r2))e-«-(*M)(l - a(r)) dr 



T2 



Now ([65ll follows from this by substituting T2 = s + r{s) and using (l55l) . 
We define (similarly to ( iSTll ) 



9 V(t x)^ 

92l^(t, X) 0<3/<x 0<y<^ 

d^V{T,x) 0<y<x 0<y<x 
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Differentiating (l50l) witfi respect to x we get 

V'(r, x) = -V'(r, xfair) - V(r, x)V"{t, x)a(r) - (1 - a(r))e~^ (68) 
V"(r, x) = -3V'(r, a^) V"(r, x)a(r) - V(r, x) V'"(r, x)a(r) + (1 - a(r))e-^' (69) 

Using this and f l67l) we obtain 



(70) 

Lemma 4. //m2(0) = Xlfcli ' "^fclO) < +oo, t/ien /or any solution of the integrated 
Burgers control problem (l43l) . ( l45l) . ( l44l) with a control function satisfying (l40l ) and for 
{t,x) G [0,t] X (0,a;] wje /iawe 

£;(t,x)xl (71) 

Proof. E{0,x) = E(0,x) x 1 follows from m2(0) < +oo. For t > we use the formula 
(COD to show that < |:E(r, 6,x(t(r))) < 3. Since < e-«'--(*M)(l - a(r)) < 1 by gTj) 
we only need to show 

< ^(3V"(x) + V'(x)) - 3^ + ^ < 3^ < 3 (72) 

The lower bound follows from 3V"{x) + V'{x) = E^i(3fc^ - k)vke-''^ > 0. 
The upper bound follows from Schwarz's inequality: 

Integrating ((TOl), using < |:E(r, ^t,^.(t(r))) < 3, ((531), dSSl) and the last inequality 
in fllOl) we obtain 

^(0, ^,.(0)) < ^(t, a;) < E{0, et,.(0)) + 3(t/2 + 1). 

Next we observe that x < ^t,^(0) < x + 1 by (l66| and -1 < \/(t, x) < 0. 
The last two bounds yield for (t, x) G [0, t] x (0, x] 



< E^{0,x + t) < E{t, x) < E*{0, x + t)+ 3(t/2 + 1) <oo. 



□ 



Lemma 5. // m2(0) < +oo, then for any solution of the integrated Burgers control 
problem (Il3l) . fl^H) . fllSll mt/i a control function satisfying (Il0|) t/iere «s a constant C* 
which depends only on the initial conditions and T such that for Tg^i < ti < t2 < T we 
have 

e{t2)-e{ti)<c*-{t2-ti) (73) 
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Proof. 6{t) = —V{t,0^). Since V(t,x) arises from (IITI) . we assume —1 < V{t,x) < 0, 
V'{t,x) < for all x > 0. 

Let us pick an arbitrary x > 0. Let C be a constant such that E(t,x) < C for 
{t,x) e [0,T] X (0,x]. 

First we are going to show that 

VO<t<T,0<a;<^ V'V{t,x) := V'{t,x)V{t,x) < C* := max{l,2C} (74) 

Note that we cannot use ( l58l ) here since that bound uses V{t,0) = 0. But 1^(0,0) = 
holds, thus ((Tl) holds for t = 0. From ([50]) and (EHD we get 



^(V'V(...)) = 

(-2V(r, x) V'(r, x)^ - V(r, x)2V"(r, x)) a(r) + (V'(r, x) - V(r, x)) e^^Xl " "(^)) < 



V(r, x)V'(r, x)M 2 - ;^V'V(r, x) ) a(r) + (V'(r, x) - V(r, x)) e-"(l - a(r)) 



From ( I5T1) we get 



V'V(r, x) > 1 ^ V'(r, x) < , < -1 < V(r, x) 

V(r,x) 

Thus by ((47]) we get 

V'V(r,x)>l =^ (V'(r,x) - V(r,x))e-^'(1-"W) < 
V'V(r, x) > 2C ^ -V(r, x)V'(r, x)^ (^2 - ^V'V(r, x) j a(r) < 

V'V(r, x) > C* =^ 4- (V'V(r, x)) < 

ar 

From V'V(0, x) < C* and the last differential inequality it easily follows by a "forbidden 
region'-argument that V'V(r, x) < C* for all < x < x and < r < T + r(T). This 
and (I53l) implies ([74]). 

By gHD and (Ell) we have 



V{ti,x) -V{t2,x) < I V{s,x)V'{s,x)ds <C* ■{t2-ti) 



for every < x < x. Letting x — 0+ implies the claim of the Lemma. □ 

3.3 No giant component in the limit 

The aim of this subsection is to prove the following proposition: 
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Proposition 2. If n^^ <^ \{n) ^ 1 and m2(0) < +00 holds for v(0) on the right-hand 
side of (l32l) then any weak limit point P of the sequence of probability measures is 
concentrated on the set of conservative forest fire evolutions: 

00 

F{Y,Vkit)^l)=l (75) 

k=l 

We are going to prove Proposition [2] by contradiction: in Lemma [6] we show that if 
9{-) ^ in the limit, then there is a positive time interval such that 9(t) has a positive 
lower bound, and that this implies that even in the convergent sequence of finite- volume 
models, a lot of mass is contained in arbitrarily big components on this interval. Than 
in subsequent Lemmas we prove that these big components indeed burn, which produces 
such a big increase in the value of the burnt mass r(-) that is in contradiction with 
E(r(T)) < 2 + E(g(T)) < 2 + T. 

By Proposition [1] the random FFE obtained as a weak limit point is almost determin- 
istic: ([371) holds with a possibly random control function r(-). Also, by (133]) we P-almost 
surely have q{t) < t from which (i40ll follows. Thus (iTTI l and (i73ll hold P-almost surely 
for the random flow obtained as a weak limit point with a deterministic constant C*. 

Lemma 6. // P„ ^ P where P does not satisfy (l75l ) on [0,T], then there exist ei, 62, 
£3 > and a deterministic t* G such that for every K < +00, every m < +00 and 

every sequence 

t* ~ El < ai < I3i < a2 < (32 < ■ ■ ■ < am < I3m < t* 
there exists an tiq < +00 such that for every n> uq and 1 < i < m we have 

P„ ( max 1 - V Vn,k{t) > £2] > £3. (76) 

Proof. First we prove that if P does not satisfy ( l75ll then there exist €1,62, €3 > and 
ei<t* <T such that 

P( M 9it)>E2)>es. (77) 
Since (l75l ) is violated, we have P( 

supo<t<T^(^) > > £ for some e > 0. 
Let L := ^nd ti := ^ for 1 < i < L where C* is the constant in ((731). Since 

^(0) = we have 

L 

e ■ 



{ sup e{t) > e] 'zyj {e{t,) > -] 



0<t<T . , 

— — 1=1 



almost surely with respect to P. Thus P(6'(t*) > |) > f for some t* G {ti, . . .t/,}. Using 
fl73D again (ETD follows with ei := €2 := |, £3 = f • 
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Now given K and the intervals [ctj, 1 < i < m we define the continuous functionals 
/, :P[0,r]^Rby 

/. (v(0), q(-), r(-)) := T (l - v,{t))dt 

where Vkif) is defined by ( fTOll . Thus for all i 

H, := {(v(0), q(-), r(-)) G I?[0, T] : /, (v(0), q(-), r(-)) > £2} 

is an open subset of PpjT] with respect to the topology of Definition [TJ Thus by the 
definition of weak convergence of probability measures we have 



lim P„(i7,) > ¥{Hi) > P inf e{t) > s^] > 

n^oo \t* —e\<t<t* I 

from which the claim of the lemma easily follows. □ 

Lemma 7. // <ti A(n) then for every 62 > there is a €4 > such that for every 
i > there is a K and an ni such that for all n>ni 1 — Yl!k=i ^n,A:(0) > £2 implies 

E„ (r„(t)) > £4 (78) 

The proof of Lemma [7] will follow as a consequence of the Lemmas [8] and [9l 

Proof of Proposition [B We are going to show that if there is a sequence P„ such that the 
weak limit point P violates ( l75i l then for some n we have 

E„(r,(T)) >T + 2 (79) 

which is in contradiction with (l35l) and (l26l) . In fact, T + 2 could be replaced with any 
finite constant in ( l79l ). but T + 2 is big enough to have a contradiction. 

We define 5i, 52, £^3 > and t* using Lemma El Next, we define £4 using this 62 and 
Lemma [71 Given these, we choose t be so small that 

We choose K and rii big enough so that (l78l ) holds. Further on, we fix the intervals 
[ttj, 1 < i < m = so that a^+i — (3i> i holds for all i and also T — [3m> i holds. 
We choose tiq such that (TfGl) holds and let n := max{no, ni}. 
Finally, we define the stopping times ri, r2, . . . , by 

K-l 

Ti := (3i A min{t : t > and 1 - ^ 'yn,fc(i) > £2}- 

fc=i 

We have n + t* < jdi + t* < a^+i < Tj+i. 

Using the strong Markov property, (TfSl) and ([76l) . the inequality (I79|l follows: 



E(r„(r)) > ^E( + t) - r„(ri) I Ti < /3i )P( < /5j ) > me4,e3. 
1=1 

□ 
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Lemma [7] stated that if initially a lot of mass is contained in big components, then 
in a short time a lot of mass burns. We prove this statement in two steps: in Lemma [8] 
we prove that if we start with a lot of mass contained in big components, then in a short 
time either a lot of this mass is burnt or the big components coagulate, so a lot of mass 
is contained in components of size -n}/^ (the same proof works if we replace the exponent 
a = 1/3 by any < a < 1/2). Then in Lemma [9] we prove that if we start with a lot of 
components of size n^/^ then in a short time a lot of mass burns. 

We will make use of the following generating function estimates in the proof of Lemma 
El If V{x) is defined as in (liB and if v G Vi then for e < | 

K-l 

l-Y,vu>e =^ V{l/K)<{e-^ -l)e (80) 
fc=i 

eK/2 

V{l/K)<-e =^ l-^^;fc>£/4. (81) 

k=l 

Lemma 8. There are constants Ci < +oo, C2 > 0, C3 > such that if 

K~l 

1 - ^"''^(0) ^ ^2 (82) 
fc=i 

for all n then 

n 

lim P ( V Vn,k (t) + r„ (i) > 6262 ) = 1 (83) 

Where i= 

Sketch proof. If we let n — cxd immediately, we get that the limiting functions vi(t) , V2{t) , . . . 
solve (l37D . ( l38ll . (I39l l with a possibly random control function r(t) = rodt). 
The n — >■ 00 limit of (l83ll is 

e{i)+r (i) > 6282 (84) 

Now we prove that if v(-) is a solution of ((37]), ((38]), ((39]) then 1 - Ef=Y^fe(0) > ^2 
implies (IMl) with Ci = 4 and €2 = \- This proof will also serve as an outline of the proof 
of Lemma [H 

In order to prove (l84l) define V{t^x) by (1411) . Thus V{t,x) solves the integrated 
Burgers control problem (l43l) . ((44]), (1451) . 

Define f/(t, x) := V{t, x) - r{t)e-\ Thus U'{t, x) = V'{t, x) + r(t)e-^ and by ((43]) we 
have U{t,x) = —V{t,x)V'{t,x). Define the characteristic curve ^(■) by 

m = v{t,m) m = ^ m 
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Let u{t) := U{t, ^{t)) - V{0, i). Thus u{0) = 0, and 



u{t) = u{t, iit)) + u'it, mm = -v{t, m)v\t, m+ 

{V'it, m) + r(t)e-5W) Vit, m) = r(t)e-«(V(t, m) < 0- (86) 



Thus u(t) < 0, moreover 



V{t, m = V{0, ^) + r(t)e-«W + u{t) < V{0, ^) + r{t), 



(87) 



m = i + r «(^) ds + r r(s)e-«Wc?s + tV(0, 1) < 1 + t ■ r(t) + tV(0, ^). (88) 



K 



K' - K 



K' 



By ([HOD we have \/(0, ;^) < -\e2. In order to prove that e{t)+r{t)> \e2 with t = 
we consider two cases: 

If r (f) > then we are done. If r (t) < define r := min{t : ^(t) = 0}. By 
we have 



1 



£2 < 



1 1 2 







Thus r < t. By dlTD we get 



-9{t) = V{t,0) = V^(r,e(r)) < -^e^ + ^^2 = -^^2 
Thus \e2 < 6{t) < ^(r) + r(r) < ^ (t) + r (t) because by (EH) the function ^(t) + r{t) is 



increasing. 



□ 



To make this proof work for Lemma [8] we have to deal with the fluctuations caused 
by randomness, combinatorial error terms and the fact that A(n) only disappears in the 
limit. 

Proof of LemmalEi Given a FFF obtained from a forest fire Markov process by fl29l) . (l30l) 
and fl3Tll. define 



Un{t,x) -.= ^2 



k=l 



k-1 



e-""" - 1 - X{n) 



i,k{0) + 2 X] ^ri,uk-i{t) - kqn,k{t) - rn,k{t) 
1=1 

By (fTOl) we have 

n 

U^{t,x)+r^{t)e~^ = ^t;„,fe(t)e-^^ - 1 - A(n) =: Vn{t,x) - \[n) =: l^„(t,x). 



fc=i 
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W'{t,x) = -J2k-Vn,k{t)e-'''' 

k>l 

1 k 
-7;dx {W{t, x) + l + X{n)f = J2oYl ^n,i{t>n,k-i{t)e 



■kx 



k>l 1=1 



w"it,x) = j2k'-vn,kity 



kx 



k>l 



k>i ^ ^ 



n2\k]-v.{t)e 



-kx 



If X{t) is a process adapted to the filtration let 

LX{t) := lim ^E{X{t + dt) - X{t) \j^t) 
Using the martingales of Proposition [U we get 



LUn{t,x) = J2 



k>l 



k 



-^Lqnj,k-i{t) - k ■ Lqn,k{t) - Lrn,k{t) 



1=1 



-kx 



k>l 



fc-1 



1=1 



- ^ ( Vn,l{t)Vn,k-l{t) 



I ■ l[2l = k] 



n 



Vn,l{t) - 



k ■ { Vn,k{t) - -1Jn,k{t) ] - (A(n) ■ k ■ Vn,k{t)) 



-kx 



- iW{t, x) + l + Xin)f - -W"{t, 2x) + 
2 n 

W'(t, x) + -W"(t, x) + X(n)W'(t, x) = 
n 

- x)Wn{t, x) + ^ x) - W:{t, 2x)) (89) 

Given the random function Wn{t,x) we define the random characteristic curve ^n{t) 
similarly to f l85ll : 

Ut) = Wn{t,Ut)), U0):=^ (90) 

This ODE is well-defined although Wn(t,x) is not continuous in t, but almost surely 
it is a step function with finitely many steps which is a sufficient condition to have 
well-posedness for the solution of (l90l) . Define u„(t) := Unit.^nif)) — W^n(0, -^). Thus 
M„(0) = and 

Un{t) = W^n(t,a(t))-W^n(0,^)-r„(t)e-«"W = K(t,^n(t))-K(0,^)-r„(t)e-«"W (91) 



The solution of ( l90l) is 



(92) 
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Putting together fl89ll and (l90ll similarly to (l86l) and using (l56ll we get 

Lur,{t) < - {W:{t,Ut)) - W:it,2Ut))) < n-' ■ Ut)-' (93) 
n 

Now Un{t) = Unit) — Lun{s)ds is a martingale and 
Lun{tf= lim \-E{{U^(t + h,Ut))-Un{t,Ut))Y\^t) < 

2 



k,i=i ^ ^ 
Define the stopping time 

Tn := min{t : .^„(t) = a = 1/3. 

In fact any < a < 1/2 would be just as good to make the right-hand side of (l93ll and 
( l94l) disappear when t < r„ and n — > oo. 

It follows from ( l94l ) and Doob's maximal inequality that 

sup |M„(t A A T)| ^ as n ^ oo 



t 



By ([93]) we have M„(t) + /g ""^ • ^^(■s) ds > Un{t) thus 

sup M„(t A r„ A T) ^ as n — > oo (95) 



t 

By (IHOD and ((82]) we have 



V;(0, ^) < (e-i - 1)62 =: -e, (96) 



Define the events A„, _B„ and the time t„ by 

/■* 1 

A„ := { sup / M„(s)rfs < — } n {M„(r„ AT) < £5/3}, 

Bn ■■= {r„(r„) < £5/3}, 



3 3 



We are going to show that that there are constants C2, C3 < +00 such that 

n 

An^{ Yl ""^'^ + ''^ ^ ^2^2} (97) 
28 



which, since (l95ll implies that lim„_,oo P ( ) = 1, gives (l83ll . 
First we show that 

AnOBnC {r„ < in}. (98) 

"in 



If we assume indirectly that An, and r„ > t„ hold then J^" Un{s)ds < so by (l92l) 
we get 

11 /"*" - 1 - 

e«(tn) + + J rn{s)e-^-^'^ds + inWniO, ^^(0)) < "^^ + " r„(r„) < 0. 

But ^n(tn) < is in contradiction with r„ > t„, thus (l98ll holds. 

Now, by (jnH) we have K(rn,n-i/3) = m„(7:„) + K(0, +r„(r„)e-""'^'. Thus by ii, 
the definition of An and and ( ISTl) we get 

An n 5„ C {M„(r„) < y } n {V;(0, -^) < -e^} n {r„(r„)e-""'' < f } ^ 

n 

Thus we have 

n 

An C {An nBn)UB'^C{ ^ Vn^'^n) > ^5/12} U {r„(r„) > £5/8} C 

fc=rii/3e5/6 

n. 

{ ^ fn,fc(T-„) + r„(r„) > Csera} 

with C3 = (1 - e"^)/6 and C2 = {1 - e-^)/12. But ZlLcseani/^ + rnit) increases 

with time, from which (l97l) follows. □ 



Lemma 9. There are constants C4 < +00, C5 > SMc/i t/zat if 

n 

for all n then with 

in ■■= Cief{n-'/Hog{n) + {n\{n))-^) (99) 

we have 

lim E(r„(y) >C5£2. (100) 



Remark. The upper bound ( 19911 is technical: on one hand it is not optimal, on the other 
hand, for the proof of Lemmal2\ we only need t „ -C 1 as n 00. 
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Proof. If f is a vertex of the graph G{n,t) let C„(f,t) denote the connected component 
of V at time t. Denote by Tb{v) the first burning time of v: 

n{v) := inf{t : \Cn{v,t+)\ < |C„(f,t_)|} 

Of course |C„(f , rb(f )+)| = 1. Define n := C^e2n^^^ and 

Unit) := {v : \Cniv,0)\ > n and n{v) > t} 

Fix a vertex v G 7Y„(0). 

n 

Wn{t) ■.= -\nn{t)\ 

n 



Znit) ■=- ^ ^{rb{w)<t} = ^^n(0) - Wnit) 



n 

u>e-H„(o) 

Thus c„(t) is an increasing process (we "freeze" c„(t) when it burns). We consider the 
right-continuous versions of the processes Cn(t) , Wn(t) , Zn{t) . 

Wn{0) > 0262/2 =: £6. 

We are going to prove that there are constants C4 < +00, C5 > such that 

limE{z^{Q) >C562 (101) 

n— >oo 

which imphes (11 001) . 

Define the stopping times 

:= mi{t : w„(t) < £6/2} 
inf{t : c„(t) > £6/4} 



^9 



Since v G 7in(0) we have 



Cn[t) > C„(0) = > 



n n 

If Cn{v,t) is connected to a vertex in Tinif) by a new edge at time t then 

c„(t+) - c„(t_) > -, log(c„(t+)) - log(c„(t_)) > log 1 + — - > ' 



n \ ncn{t-) J ncn{t-) 

loef2)n 1 

\Cn{v,t)\ {\nn{t)\ - \Cn{v,t)\) l{i > log(2)n ■ (m;, (t) - Cn{t)) > 

log(2)n|l|,<,} = n^/^i^ ■ C2 ■ C3 ■ (£2)' ■ =: n'/'e,l{t<r} 
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Thus log(c„(t)) — St n^^'^{t /\t) is a submartingale. Using the optional sampling theorem 
we get 



n 



-£7-ni/'E(r) >E(log(c„(r))) -£7-^'/'E(r) > log(c„,(0)) > - log( 
By Markov's inequality we obtain that for some constant C < +00 

P(r<Cn-i/VlogW) >^ 
If < T6('i^)Ar^„, then C„(f,rg) > ^n, so E( rfe(t>) — ) < (nA(n))~"'^^, which implies 

P(r^ ATft < Cn-^/^£2^1og(n) + C"(nA(n))-V2"^) > ^. 

for some constant C . We define i of (l99l ) with C4 := max{C, C'}. Using the linearity of 
expectation we get 

E(2;„(t)) =E(^ lKM<ti) >£6P(T-5(i;)<t). 

The inequality l{T„<f}^ < follows from the definition of Tyj. 

7 < P(r^ ATfe <t) < P(r^ <t) +P(r, < t) < E( ) - + E( ) - 
From this ffTnTD follows. □ 

4 The critical equation 
4.1 Elementary properties 

Existence to the solutions of (i37l ). ( l39l ) with initial condition satisfying 777.2(0) < +00 and 
boundary condition 

00 

Y,^k{t) = l (102) 

k=l 

follows as corollary to Propositions [U and EJ indeed for any initial condition vq G Vi 
we can prepare a sequence of initial conditions of the random graph problem such that 
(l32l) holds as 77 — > cxD (we do not need to assume convergence of 7r7„^2(0) to 7772(0)). If 
77^^ <^ A (77) <^ 1 then any weak limit of the probability measures is concentrated on 
a subset of FFFs which generate a FFE satisfying f l37ll , p02p . 

Moreover it is easily seen that p02p implies that r(-) must be continuous, and for 
fc > 2, the functions t 1— > Vk(t) solving ( l37l l are differentiable. Thus v(-) solves (fT2l ). f |T3ll . 

Note that assuming that v(-) G Svo[0,T] is a solution of (fT2l ) . (fT3ll one can deduce 
only from these equations that (1371) holds with a control function r(-) satisfying (llOll : 
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one has to define a FFF using (l33l) and qk,oo{-) = 0: plugging 9{t) =0 into (!2T|1 we can 
see that the function r(-) is increasing. 

Taking the generating function of a solution of (l37l) . (l39l) . (1102^ with initial condition 
satisfying m2(0) < +cx3 we get a solution of ( l43ll . (i45l) satisfying the boundary condition 
l^(t,0) = 0. 

In this case the increasing function t r{t) is absolutely continuous with respect to 
Lebesgue measure: its Radon-Nykodim derivative r(t) = ip{t) is a.e. bounded in compact 
domains: 

Taking the limit x — in (i43l l and using (iTTI ). ( l58l l (which holds because V{t, 0) = 0) 
we find ^ 

r(t2) - = lim ^ / x)ds <C-{t2- h). (103) 

x^o 2 

Thus in the sequel we assume given a solution of the critical Burgers control problem 

V{t^ x) = -V'{t, x)V{t, x) + e-Xt), (104) 

V{t,0) = (105) 

V{0,x) = Voix) (106) 

where ip(t) is nonnegative and bounded on [0,T], and V(t,x) is of the form ( l4Tll . 

Lemma 10. For an?/ solution of 01040 . pOGf l. pOSp wi^/i V'(0) < +cx) and /or any 
t > T'gei (^see ^) we have V'{t,0) := lim^^o ^'(^, a;) = — oo. 

Proof. We actually prove that for any t < oo, x < oo there exists a constant C = 
C(t,x) > such that for any (t, x) G [Tge^t] x (0,x], -V'{t,x) > C/y/x. 

One can prove the upper bound of (iGOll for all V{x) satisfying V^(0) = without the 
assumption (l59ll (the same proof works). 

From ([7T]) and the upper bound of (l60l) it follows that there exists a constant C < oo 
such that for (t, x) G [Tgei,t| x (0,x] 

^(t, x)-^ < C, -V{t, x) < Cx^l'^. 

Differentiating with respect to x in (11041) we get 

I {-V\t, X)) = V'{t, xf + \/(t, x)V"{t, x) + e^Xt) = 
^'(t, xf . [l - ^^^^^1^) + e->(t) > V'(t, x)^ (l - C^x'/' . i-V'it, X))) (107) 

There exists a < C such that for x G (0, x] we have 

-\/'(Tgei,x)>a/v^ (108) 
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by (IBTl) and (17T1) . since V"'(Tgei,0) = — oo <^==^ mi(Tgei) = +00 follows from the fact 
that for t < Tgei the solutions of ([6]) and (fT2l ) + (fT3ll coincide, and it is well-known from 
the theory of the Smoluchowski coagulation equations that we have ([9]) for the solution 

of (El). 

From the differential inequality (11071) it follows that 

- V'{t, x) < Ix-^/^ =^ ^ {-V'{t, x)) > (109) 

Let C := min{C, C^^}. For (t,x) G [Tgei,t] x (0,x] the inequality 

-V'{t,x) > C/v^. 
follows from (llOSp and pOOp by a "forbidden region"- argument. 



□ 



Summarizing: from Lemmas [HI (H [lO] and (11031) it follows 
Lemma 11. For {t,x) G [Tgei,t\ x (0,x] 



-V{t,x) 




-V'{t,x) 


X X-'/' 


V"it,x) 


X X-'/' 


V{t,x)V'{t,x) 


xl, 




X 1. 



(110) 

(111) 

(112) 
(113) 
(114) 



4.2 Bounds on E' 



In this subsection we assume given a solution of ( 1104p . pOSp . (11061) satisfying |V^"'(0, 0)| < 
+00. All of the results of the previous subsection are valid for V{t,x). 

Lemma 12. 

E'iTg,t,x) = Oix-'/^) (115) 

Proof. We consider the function X{t,u) defined for every t as in the proof of Lemma [3l 
X"'{0,u) = 0(1) for u G [0,u] by mi(0) > and m3(0) < +00. For t < Tgei we have 
ip(t) = thus V(t, x) satisfies the Burgers equation 

V{t,x) + V{t,x)V'{t,x) =0 

from which 

X{t,u) = X{0,u)-tu 
follows. Differentiating ( l63l l with respect to x we get 

E'(Tgei,a;) = E(Tgei,x)=^X'"(0,-\/(Tgei,x))\/'(Tgei,x). 
Now (fTT5H follows from ((711) and (EI]). □ 
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From now on, we consider the solution of (11041) . (11051) . (I106p for t > Tgei, that is we 
assume that Tgei = 0. 

Since the function r(t) is continuous we get that t(r) defined by (l46l) is the inverse 
function of t + r{t) which by (1481) implies V(r, a;) = V(t{T),x). Integrating (iTOl) and 
using (l53l ). ( I55l) we get for < ti < t2 < C)0 



(116) 
(IIT) 

Lemma 13. The function (t,x) i-^ E{t,x) is continuous on the domain (t,x) G [0,t] x 
[0,x], and 

(fi{t) = lim V\t, x)V{t, x) = E{t, 0). (118) 

Proof. From (11141) and (l65l) it follows that the characteristic curves ^t,x{s) are jointly 
continuous in the variables {(t,x,s) : < s < t, < a;}. And hence, further on, 
from (11161) and ((72]), by dominated convergence it follows that (t,x) i— > E(t,x) is jointly 
continuous in {{t,x) : < t, < x}. Further, from (l62l ) it follows that 

limV{t,x)V'{t,x) = lim E{t,x) =: E{t,0) 

x—^O X— »0 

Hence, pi8p follows from (11031) again by dominated convergence. □ 
Lemma 14. 

(i) The function x i-^ E{t,x) is Hdlder-1/2 at x ^ 0: 

Eit,x) = ip{t){l + 0{x'/^)). (119) 

(a) The function t t-^ ip(t) is Lipschitz continuous: there exists a constant C < oo 
(which depends only on the initial conditions p06p and the choice oft such that for 
any ti,t2 E [0,t\ 

\^{ti) - y^{t2)\ < C\ti - t2\. (120) 

Proof (i) We prove \E'(t, x) \ = C(x^^/^). In this order we shall use the following a priori 
estimates 

6,.(s)x (xi/2 + (t-s))' (121) 
Cxi') ■■= dx^tAs) = O {{x"' + {t- s))x-'/') . (122) 
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Indeed: ffml) follows from ((65]), (fTTnl) and (fTTil) . and we get (fT22ll from (fTTTll and from 
the fact that characteristics do not intersect (thus < Ct,xi^)) by differentiating (l65l) 
w.r.t. x: 

o<CA-'^)<i-v'it,x){t-s) 

The a priori bound 

\E'{t,x)\ = 0{x~^). (123) 

follows from 

E'{t, x) = -3V'it, xf + E{t, x) = 0((a:-^/^)^) + 0(x-^) 

by dm]), dnl) and 

-^v'"(t,x) < r r'"(y)rfi/ < v"(^) = o{x-'/') 



using both the upper and lower bounds of (11121) . 
Differentiating with respect to x in (11171) yields 



E'(t,x) = i?'(0,6,x(0))el.(0)+ (124) 
/•* . E'( 



,^'(s,6,x(s)) , 3i?(s,6,x(s))r'(s,6,.(s)) 



(s,6,.(s)) V'(s,6,.(s))2 

+V'(.,6,.(.)) v'{s,^us)r^^'''^^ 

Next using (I123P bound we estimate the expression of E'{t,x) given in (I124p . Using 
dni), (fTTTll . (ITtId . (fTT5D . (fT2TD . and (fT22D we conclude that if (fT^H]) holds then actually 

\E'{t,x)\=0{x-^^^). (125) 

The dominating order is given by the first term (outside the integral) and the first two 
terms under the integral on the right hand side of p24p . 
Finally, (fTTOll follows from ( fTTSD and (fT25D . 

(ii) In order to prove (11201) we note that from (11161) and (IllSp it follows that for 

0<ti<t2<t 

^{h) - ^{t2) = E{h, 0) - E{ti,^t„o{ti)) 



s 



Hence, by (fTTOD . (fT2TD and (CSD we obtain directly (fT2nll . □ 
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Summarizing again, from Lemmas [3l (H [IHl [13] and [TH it follows 

Proposition 3. For a solution of (11041) . f llOGp . (llOSp with initial condition satisfying 
Tgei = 0, jZI]) and ffUBl) and /or {t,x) E [0,I\ x (0,x] 

-F(t, x) = ^/Mt)x'/^1 + Oix'/^)) , (126) 

(127) 



-nt,x) = y^x-V2(l + 0(xV2)). 



nt,^) = V4^^-'/'(l + (^(x^/')), (128) 

V{t,x)V'{t,x) = ^{t){l + 0{x'/^)). (129) 

V{t,x) = 0{x^^^), (130) 

V'(t,x) = 0(x-^/2)^ ^131^ 

(/.(t)xl, |</^(ti)-</^(t2)| <C|ti-t2|. (132) 

In order to prove f |T4l) we need Example (c) of Theorem 4. of chapter XIII. 5 of 0. 
With our notations each of the relations 



-V{t,x)r^x'-'/'^/2^) and ^^,(t) ~ ^fc^/^-i^/^ 
implies the other. 



l=k ^^2' 



4.3 Uniqueness 

We are going to prove Theorem [TJ by proving the uniqueness of (I104p . (I106p . (11051) . 

Proof of TheoremUi Assume that V(t,x) and U{t,x) are two solutions of the critical 
Burgers control problem with the same initial conditions and with the control functions 
ip{t) and 'ipif)^ respectively. Denote 



*(*):= (134) 
Then, it is easily seen that that (given S{t,x)) W{t^x), 5{t) will solve the linear control 



36 



problem 

W{t, x) + {S{t, x)W{t, x))' = e-^5(t), (135) 

W{0,x)=0, (136) 

W{t,0) = 0. (137) 

We assume S(t, x) and p(t) given, with the regularity properties inherited from Proposi- 
tion [3j 



-S{t, x) = ^/Mt)x'/^ (1 + Oix'/^)) , (138) 
-S'it,x) = ^x-'/'{l + 0{x'/')), (139) 



(140) 



S"it,x) = ^x--'/'{l + Oix'/')), 

S{t, x)S'{t, x) = p(t) (1 + 0(x^/2)) . (141) 

S{t,x) = 0{x^/^), (142) 

S'{t,x) = 0{x-^/^), (143) 

p(t)xl, |p(ti)-p(t2)| <C|ti-t2|. (144) 

We will prove that under these conditions, the unique solution of the problem (I135p . 
ffT36D . ffnTI) is W{t,x) = 0, 5{t) = 0. 

First we define the characteristics of the equation 01351) : these are the curves [0,t] 3 
s ^ Ctis) defined by the ODE 

Us) = S{s, Ct{s)), Ct{t) = 0, Ct{s) > for s < t. (145) 

Next we define the functions [0, t] 9 s i— > Pt{s) 

f3t{s) ■.= S'{s,Ct{s)). 

The functions [0,t] 9 s Ct{s) and [0, t] 9 s ^— Pt{s) are directly determined by S{t,x) 
and from 01380 . (11391) . OMOp and 0144^ inherit the following regularity properties to be 
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used later: 

as)= ^{t-sr{i+o{t-s)), (146) 

Q{s) = -p{t){t-s){l + 0{t-s)), (147) 
Ctis)= p(t){l + Oit-s)), (148) 
f3,{s) = -{t-sr'{l + 0{t-s)), (149) 
$tis) = -it-sr^{l + 0{t-s)). (150) 
We define [0, t] 9 s i— > r]t(s) as 

r]tis):=Wis,as)), 

with W{t,x) given in (fT34ll being solution of ( fT35D . (fT36ll . (fT37ll . Then, for any t > 0, 
5(s),?7t(s), s G [0,t] solves the ODE (boundary value) control problem 

litis) + PtisMs) = e-^'(^)5(s), rit{0) = = r]t{t) (151) 

We will prove that this implies 6{t) = 0. Hence it follows that W{t,x) = 0. 
On the domain {(t, s) : < s < t < oo} we define the integral kernel 

t - s 

/C(t, s) := exp { / Ptiu)du- Ct{s)} = ——C{t,s), 
Jo 't 

defined on the same domain {(t, s) : < s < t < oo}, where 

C{t, s) ■= exp { / (A(n) + {t- u)'')du - 0(s)}. 



The ODE control problem p5ip is equivalent to 

/C(t,s)5(s)ds = 0. (152) 
It is handy to introduce the function 

7(t) := /" S{s){t - s)ds. 



Then, after two integrations by parts the identity (I152p is transformed into the eigenvalue 
problem 

-t 

£(t,s)7(s)rfs = 7(t), (153) 
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, 2dsC{t,s)-{t-s)dlC{t,s) 



where 

Z{t,s):={dMt.t))-'dlK:{t,s, ^^^^^ 

Using the regularity properties (fT46ll . f fT47l) . f fTiSD . f fT49D . f fTBOD it follows that 



sup 

0<s<t<i 



/C(t,s) < cx). (154) 



From (I153P and (I154p . by a Gronwall argument we get 7(t) = and hence 5{t) = = 
W{t^x), which proves uniqueness of the solution of (I104p . (I106p . (llOSp . 

□ 



References 

[1] D. Aldous: Deterministic and stochastic models for coalescence (aggregation and 
coagulation): a review of the mean-field theory for probabilists. Bernoulli, 5: 3-48 
(1999) 

[2] J. van den Berg, R. Brouwer: Self-destructive percolation. Random Structures and 
Algorithms, 24: 480-501 (2004) 

[3] J. van den Berg, R. Brouwer: Self-organized forest-fires near the critical time. Com- 
munications in Mathematical Physics, 267: 265-277 (2006) 

[4] J. van den Berg, A. Jarai: On the asymptotic density in a one-dimensional selfor- 
ganized critical forest-fire model. Communications in Mathematical Physics, 253: 
633-644 (2005) 

[5] J. van den Berg, B. Toth: A signal-recovery system: asymptotic properties, and 
construction of an infinite- volume limit. Stochastic Processes and their Applications, 
96: 177-190 (2001) 

[6] B. Bollobas: Random Graphs. Cambridge University Press, 2001 

[7] William Feller. An introduction to probability theory and its applications. Vol. II. 
Second edition. John Wiley & Sons Inc., New York, 1971. 

[8] R. Brouwer: Percolation, forest-fires and monomer-dimers (or the hunt for self- 
organised criticality) . PhD thesis, VU Amsterdam, 2005. 

[9] E. Buffet. J.V. Pule: On Lushnikov's model of gelation. Journal of Statistical 
Physics, 58: 1041-1058 1990 

[10] E. Buffet. J.V. Pule: Polymers and random graphs. Journal of Statistical Physics, 
64: 87-110 1991 



39 



[11] B. Drossel, F. Schwabl: Self-organized critical forest fire model. Physical Review 
Letters, 69: 1629-1632 (1992) 

[12] M. Duerre: Existence of multi-dimensional infinite volume self-organized critical 
forest-fire models. Electronic Journal of Probability, 11: 513-539 (2006) 

[13] M. Duerre: Uniqueness of multi-dimensional infinite-volume self-organized critical 
forest fire models. Electronic Communications in Probability, 11: 304-315 (2006) 

[14] P. Erdos, A. Renyi: On random graphs I. Publicationes Mathematicae Debrecen, 6: 
290-297 (1959) 

[15] S. Janson, T. Luczak, A. Rucinski: Random Graphs. John Wiley and Sons, NY, 
2000 

[16] K. Schenk, B. Drossel, F. Schwabl: Self-organized critical forest-fire model on large 
scales. Physical Review E, 65: 026135, (2002) 



40 



